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Abstract 

Using over 40 million muons collected since 1989 by the MACRO detector we have searched for a depletion 
of muons coming from the direction of the Moon due to primary cosmic rays striking the Moon. We observe 
this Moon shadow in the expected position with a statistical significance of more than 5 standard deviations. 
We have analyzed the same data for an analogous Sun shadow, and have found a signal with a significance of 
about 4 standard deviations. The Sun shadow is displaced from the Sun's position by about 0.6 degrees North 
in ecliptic coordinates. This displacement is compatible with a deflection of primary cosmic rays due to the 
Interplanetary Magnetic Field in the 10-20 TeV primary energy range which is relevant to the underground 
muons observed by MACRO. 

x. 

u ; 1 Introduction: 

MACRO was primarily designed to search for magnetic monopoles and rare particles in the cosmic rays, 
including high energy atmospheric neutrinos and muons from cosmic point sources (Ahlen et al., 1993). An 
important goal of the MACRO detector therefore is the recognition of point sources of high energy cosmic 
t— i ' rays (E> 10 TeV) by looking for an excess of underground muons from a particular direction of the sky and 
above a nearly isotropic background of cosmic rays. As previously suggested (Clark, 1957) the shadow in the 
\ high energy cosmic ray flux from the directions of the Moon and the Sun could be observed and employed to 
measure the angular resolution and pointing accuracy of the detector. Clear evidence of this effect has been 
reported by many EAS experiments (Aglietta et al., 1991, Alexandreas et al., 1991, Amenomori et al., 1993, 
! Borione etal., 1994). 

Q\ | MACRO presented the first deep underground evidence of the Moon shadow effect (Ambrosio et al. 1999) 

and now presents evidence of the Sun shadow due to its excellent angular resolution (ps 1°) and its large col- 
lecting area and high statistics when compared to other deep underground detectors. Another very interesting 
point is that almost all cosmic ray particles are charged and deflected by Geomagnetic (GMF) and Interplan- 
etary (IMF) magnetic fields. The shadows of the Moon and the Sun show these effects by means of displaced 
obscurations, whose analyses might give new information about these fields as well as about the cosmic-ray 
. i-h ! energy and charge distribution. 

For this analysis we use the same technique used by Alexandreas et al., 1991, Amenomori et al., 1993 
and Borione et al., 1994. In this work we extend the sample used for the previous Moon shadowing analysis 
(Ambrosio et al. 1999) and present the results for the Sun shadowing effect. 

2 Event Selection: 

The muon sample used in this study includes all events recorded between February 1989 and December 
1998. About 45 x 10 6 events have been collected over 2.6 x 10 3 live days. The criteria used to select events 
for this analysis are the same as in Ambrosio 1999, defined to optimize the quality of reconstructed tracks. 
The selected events are consequently those which most accurately point back to their origin on the celestial 
sphere. These selection criteria reduce the sample size to 38.2 x 10 6 muons. 

The topocentric positions of the Moon and the Sun are computed at the arrival time of each event in the 
sample, using the database of ephemerides available from the Jet Propulsion Laboratory, JPLEPH (Standish 
E.M. et al., 1995) and include a correction for the parallax due to MACRO'S instantaneous position on the 
earth. The muon events in two 10° wide windows centered on the Moon and the Sun are selected for further 
analysis. 



There are 3.17 x 10 5 events in the Moon window that pass all cuts. For the Sun shadow analysis there are 
3.08 x 10 5 events collected that pass all cuts. Twenty five background samples are generated for each run used 
in the analysis. These backgrounds are constructed by coupling the direction of each muon in the run with 
the times of 25 randomly selected muons from the same run. The 25 background samples are then processed 
using the same procedure as the muon data sample. 

3 Shadow of the Moon: 

The search for the Moon shadow in a direction- 
independent way and the estimate of its signif- 
icance is possible with the maximum likelihood 
method of COS-B, a technique first described in 
detail by Cash W., 1979. This method is based on 
a priori knowledge of the point spread function of 
the MACRO detector (PSF). We have determined 
the PSF using the observed space angle distribu- 
tion of double muons divided by y/2. The PSF, 
as determined from more than 10 6 muon pairs, 
looks strongly peaked and has a non-Gaussian 
behavior (Ambrosio et al. 1999). To find the 
most likely position of the Moon, we compare 
the two-dimensional distribution of muons in the 
window centered on the Moon with the expected 
background events in the same window. In this 
analysis, each muon event is first sorted into a 
grid of equal solid angle bins (AO = 0.125° x 
0.125° = 1.6 x 10" 2 deg 2 ). The shadowing 
source of strength Sm at a fixed position (x s ,y s ) 
that best fits the data is then found by minimiz- 
ing 
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Figure 1: The two-dimensional distribution of A in bins of equal solid 
angle in the Moon window. The axes are offsets from the Moon center. 
The fiducial position of the Moon, at position (0,0), is marked by a +; 
a circle corresponding to the average lunar radius, 0.26°, is centered at 
this position. The A grey scale is given at the right margin of the figure. 
The maximum of this distribution has A = 30.0. 
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where the sum is over all bins in the window. Here jVj is the number of events observed in each bin i, Nf x is 
the expected number of events in bin i, and nun is the number of bins in the grid. This expression assumes 
that a Poissonian process is responsible for the events seen in each bin. The expected number of events in bin 
i is given by Nf x = N^ kd — Sm ■ V(x s — Xi, y s — yi), where N^ kd is the average number of background 
events at the position (xj, yi), and Sm • V(x s — Xi, y s — yi) is the number of events removed from bin i by 
the shadow of the Moon. Here V(x s — Xi,y s — yi) is the MPSF, modified for the finite size of the Moon's 
disk, computed at the point (xi,yi) when the shadowing source is at (x s ,y s ). Finally, the shadow strength 
5a/ that minimizes x 2 is computed for every grid point in the window. This minimum x 2 { x s, y«> SW) is then 
compared with x 2 (0) f° r th e null hypothesis that no shadowing source is present in the window (Sm = 0), 

X = X 2 (0)-X 2 (x s ,y s ,S M ). 

In Figure^we show the results of this analysis in a window 4.375° x 4.375° centered on the Moon. This 
window has been divided into 35 x 35 cells, each having dimensions 0.125° x 0.125°. In this window the 
number of muon events is of 14388. In this figure, A is displayed in grey scale format for every bin in the 
Moon window. 

Also shown is the fiducial position of the Moon and a circle centered at this position corresponding to the 



average lunar radius, 0.26°. The maximum A = 30.0, corresponding to 5.5a, is found at the expected Moon 
position and this provides further confirmation of MACRO'S absolute pointing. The value of the shadow 
strength obtained by the likelihood method at this position, Sm = 243^} events, agrees well with the ex- 
pected value of 215 events. We have verified the properties of the A distribution by constructing 71 other win- 
dows similar to the Moon window, each displaced from the next by 5° in right ascension. For each off-source 
window, we have followed the procedure used for the Moon window in computing the expected background. 
To avoid edge effects associated with a source near the edge of a window, we only have evaluated A for the 
central 12 x 12 bins. 

4 Angular Resolution of the MACRO Apparatus: 

If we assume that the space angle distribution of double muons is a good approximation of PSF, the angular 
resolution of the apparatus can firstly be estimated by this distribution as the angle 6 G8 % of the cone that 
contains the 68% of the events from a point-like source. Using this definition of the angular resolution, we 
obtain 6 m % = 0.8°. Since a simple Gaussian function cannot fit the PSF, we cannot compute the angular 
resolution as the a of the Gaussian function that maximizes the likelihood function for the detection of the 
Moon shadow. Therefore we define a scale parameter T that rigidly scales the PSF by the factor T, V (x s — 
Xi,y s — Ui\ J 7 ) and then repeat the likelihood analysis in the Moon window for different values of T . We 
assume that the value of T that maximizes A gives the best T for computing the angular resolution. Using 
these procedure we find that the angular resolution is 0.86°7_g26° fully consistent with the other estimation. 

5 The Sun Shadow: 

We have repeated the same analysis by 
using events in the Sun window. In eclip- 
tic coordinates (ecliptic latitude and eclip- 
tic longitude), the maximum-likelihood 
method gives a map of A, in the Sun win- 
dow, as shown in Figure |2 The most prob- 
able position of the center of the deficit is 
found at (0° ,0.625°) with a A=17.4 corre- 
sponding to a significance of 4.2 a. The 
probability to find random fluctations hav- 
ing ap(x 2 > 17.4) is less than 10~ 5 , hence 
the observed Sun's shadow is significant. 

The large displacement of the shadow 
from the apparent position of the Sun is 
likely due to the combined effect of the 
magnetic field of the Sun and the geomag- 
netic field. However a quantitative inter- 
pretation of the Sun's shadow in connec- 
tion with the effect of the solar magnetic 
field is difficult, because during almost ten 
years of MACRO data acquisition, the ef- 
fect of the IMF is averaged over its changes 
due to the variable inclination and over the 
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Figure 2: The two-dimensional distribution of A in bins of equal solid angle 
in the Sun window. The axes are offsets from the Sun center using ecliptic 
coordinates. The fiducial position of the Sun, at position (0,0), is marked by 
a +; a circle corresponding to the average solar radius, 0.26°, is centered at 
yearly variation of solar activity with a pe- this position. The A grey scale is given at the right margin of the figure. The 
riod of 22 years. Moreover the IMF has a maximum of this distribution, A = 17.4, is about 0.6° N . 
sector structure with the magnetic field direction reversed so that in some sectors the magnetic field points 
inwards and in others outwards and the neutral sheet separating the two regions is inclined to the equatorial 



plane. As the earth is within 7° of the solar equatorial plan, the polarity of the IMF will change several times 
during the Sun rotation of 27 days on its axis. 

A careful selection of the data as a function of IMF polarity configuration, as well as higher statistics, could 
be used to obtain new information about the IMF. 

6 Conclusions: 

The MACRO detector, operational since February 1989, has collected a muon sample of about 40 million 
events. Using this sample we have observed the Moon and Sun shadows in the cosmic ray sky at primary 
energies « 10 - 15 TeV. In the deficit analysis, we find an event deficit around the Moon of about 5.0a 
significance. With a maximum likelihood analysis, we confirm the detection of the Moon's shadow with a 
significance of 5.5a. Our estimate of the angular resolution is # 68 % =0.86°+Q2go. These results characterize 
MACRO as a muon telescope by confirming MACRO'S absolute pointing ability and by quantifying its angular 
resolution. This investigation shows that the MACRO detector has the capability of detecting signals from 
point-like sources by detecting secondary underground muons . The analysis of the deficit in the Sun direction 
and in ecliptic coordinates gives a Sun shadow with a significance of 4.2 a. The shadow of the Moon and 
the Sun is shifted due to the curvature of cosmic rays by GMF and IMF at the primary energy of 10-20 TeV 
relevant to MACRO. 
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